Effective potential for Polyakov loops from a center symmetric effective theory in 

three dimensions 
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We present lattice simulations of a center symmetric dimensionally reduced effective field theory 
for SU(2) Yang Mills which employ thermal Wilson lines and three-dimensional magnetic fields as 
fundamental degrees of freedom. The action is composed of a gauge invariant kinetic term, spatial 
gauge fields and a potential for the Wilson line which includes a "fuzzy" bag term to generate 
non-perturbative fluctuations. The effective potential for the Polyakov loop is extracted from the 
simulations including all modes of the loop as well as for cooled configuration where the hard modes 
have been averaged out. The former is found to exhibit a non-analytic contribution while the latter 
can be described by a mean-field like ansatz with quadratic and quartic terms, plus a Vandermonde 
potential which depends upon the location within the phase diagram. 

PACS numbers: 12.38.-t, 12.38.Gc, 12.38.Mh 
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I. INTRODUCTION 

QCD at temperatures T ~ 200 MeV-1 GeV exhibits 
only partial deconfinement and, moreover, perturbation 
theory fails to reproduce some thermodynamic quantities 
such as the equation of state [1]. Several authors [2-5] 
have suggested that in this regime a more appropriate 
effective description is not in terms of quasi-particles but 
in terms of the thermal Wilson line 
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and the spatial components of the gauge field. The op- 
erator Aq from eq. (1) is defined on four dimensional eu- 
clidean space-time while L(x) is a matrix-valued field in 
space. Zr denotes a renormalization constant. In such 
a framework, the deconfined phase is not a free gas, but 
rather a condensate of spin-like operators 



f 



O) = ^trL(x) 



(2) 



called Polyakov loops. The volume averaged expectation 
value of this operator is an order parameter for the de- 
confining phase transition in the limit of infinitely mas- 
sive quarks [6]. In contrast to ferromagnetism, the high 
temperature phase is the ordered phase where a global 
symmetry is spontaneously broken. 

We perform lattice simulations of an effective theory 
in three dimensions defined in the continuum by [4] 
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All fields in (3) are functions of x only. Gij is the mag- 
netic field strength. The second term is the contribution 
from electric fields since in the three-dimensional theory, 
for arbitrary Aq, E is given by [4] 
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The potential ~ —J2 n >i |trL"| 2 /n 4 for the Wilson line 
L is obtained by computing the one-loop fluctuation de- 
terminant in a constant background Ao (or L) field [7]. 
It is evidently minimized by the perturbative vacuum 
(L) = 11 (times a phase), for any T. To generate a phase 
transition in infinite volume, refs. [3, 4] suggested to 
add non-perturbative contributions such as B/T 2 |trL| 2 , 
with Bf a "fuzzy" bag constant (see, also, refs. [8]). 
The "fuzzy bag" dominates at sufficiently low temper- 
ature and induces a transition to a confined phase with 
(trL) = 0. 

The three-dimensional effective theory (3) should be 
valid only over distance scales larger than l/T. A related 
3D theory has been formulated in refs. [5] (also see [9]). 
These effective theories respect the global Z(N) cen- 
ter symmetry of the four-dimensional Euclidean SU(N) 
Yang-Mills theory. This allows for non-perturbatively 
large fluctuations of the Wilson line and it is interest- 
ing to investigate whether such an extension of high- 
temperature perturbation theory is sufficient to describe 
the properties of hot Yang-Mills even close to the tem- 
perature for deconfinement. 

We note that there have been various attempts 
at extracting through numerical simulations a three- 
dimensional effective action which reproduces the long- 
range properties of the underlying four-dimensional 
Yang-Mills theory [10]. Here, our approach is different 
(mostly because we do not aim at matching the couplings 



2 



of the 3d theory yet). We shall focus on extracting an 
effective potential and analysing its structure from the 
3d theory itself. 

II. LATTICE ACTION 

Our present simulations have been performed for gauge 
group SU(2). The structure of this group is simpler than 
that of SU(3), and thus allows for much higher numeri- 
cal precision but exhibits the qualitative features which 
we are interested in, namely a deconfining phase tran- 
sition and non-perturbative fluctuations between center- 
symmetric states. The lattice action is of the form 

S = /?]T(l-iRctrU n ) 

□ 

-^£ tr (LiU«Ltut. + h.c.) -m 2 J2 It^l 2 • 

(ij) 1 

(5) 

The first term is the standard Wilson action for the mag- 
netic fields in three dimensions; the sum runs over all 
spatial plaquettcs. We have checked that our implemen- 
tation reproduces the plaquette expectation values pub- 
lished in ref. [11]. The second term is a kinetic term 
for the Wilson line corresponding to the electric fields 
in three dimensions. Here, the sum runs over all links 
connecting nearest neighbor sites and the gauge links 
Uij ensure gauge invariance. The third term is a mass 
term for the trace of the Wilson line which combines 
the n = 1 contribution to the one-loop potential with 
the non-perturbative "fuzzy" bag contribution. Contri- 
butions from larger n have been dropped. We have previ- 
ously performed simulations of a simplified version of (5) 
without magnetic fields in ref. [12]. For m 2 = 0, and 
without magnetic fields, our code reproduces free energy 
measurements from refs. [13] but differs slightly from the 
older work of ref. [14] which used smaller lattices and 
lower statistics. 

We employ the standard Metropolis algorithm to gen- 
erate a thermal ensemble of configurations. The lattice 
is updated sequentially. Trial steps are taken in phase 
space by rotating single sites or links by a small an- 
gle in a random direction. The difference of the ac- 
tion AS = Si+i — Si of the new and the old configura- 
tions is calculated and the step is accepted with prob- 
ability p = min [1, exp(— AS)]. To reduce autocorre- 
lations and accelerate thcrmalization we include over- 
relaxation sweeps where the Metropolis trial steps are 
taken deterministically in a way that approximately in- 
verts the action with respect to its minimum [15] (exact 
non-stochastic over-relaxation is not possible due to the 
non-linear term in the action). 

All sample sizes quoted here are taken to be statisti- 
cally independent measurements. We have estimated au- 
tocorrelation times using the binning method described 



in ref. [16] and by observing the thermal relaxation within 
the Monte Carlo time series of measurements. We shall 
refrain from discussing these technical matters here in 
detail. 



III. RESULTS 



A. Phase diagram 



To determine the phase diagram of the theory we mea- 
sure the expectation value of the Polyakov loop £(x) = 
itrL(a;) and the inverse correlation length = l/£ as 
functions of the couplings m 2 and (3. There is a line of 
phase transitions in the (3 — m 2 plane separating the re- 
gion where the expectation value of the Polyakov loop 
vanishes (the Z(2) symmetric phase at low /?, m 2 ) from 
the Z(2) broken phase at large (3 or m 2 (see fig. 1). The 
transition is of second order since the inverse correlation 
length 77i£ on the phase boundary extrapolates to zero 
in the infinite volume limit [12] as long as (3 is not too 
large 1 . It appears that the order of the phase transition 
changes to first order at roughly (3 > 3.0. This may sig- 
nal a breakdown of this model as an effective description 
of 4D Yang Mills in the extreme weak coupling limit (for 
a discussion of this issue see Appendix B). 

Confinement is realized in distinct ways. At small (3 
and vanishing m 2 we find (£} = 0.0 because the Wilson 
line li(x) averages to zero from random fluctuations over 
the group manifold. For large (3 and negative m 2 (cor- 
responding to the upper left region in fig. 1), we find a 
non-trivial Z(2) symmetric phase, where L(a;) = ir^ or 
rotations thereof, which implies that the trace of L van- 
ishes while the field L is non-zero 2 (see, also, ref. [12]). 



1 The second order phase transition is further confirmed by the 
potentials described in later sections 

2 In other words, the distributions of the eigenvalues Ai and A2, 
which are gauge invariant, peak about ±1. That distinct confined 
phases with different eigenvalue structure can also arise in 4d 
models of Polyakov loops coupled to gauge fields, was shown in 
refs. [17] 
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FIG. 1: The phase boundary in the (3 - m 2 plane. The Z(N)- 
broken phase corresponds to the upper-right region. 



to be confused with a Z(2) background field ~ — hi cor- 
responding to (heavy) dynamical quarks in the funda- 
mental representation. 

All measurements presented here were performed on a 
N s = 24 cubic lattice. However, using smaller lattices we 
have checked that the coefficients a, b, c from (9) do not 
change much with volume if N s > 12. N s — 24 appears 
to be a good approximation to the infinite volume limit 
where, as indicated in eq. (8), one expects the potential 
to be volume independent. 5000 independent lattice con- 
figurations were generated for each combination of (3, m 2 . 

We first consider the case without potential, m 2 = 0. 
The phase transition occurs along a vertical line in the 
phase diagram at (3c ~ 0.9. Below the phase transi- 
tion point we find that the potential defined via eq. (7) 
coincides with the Vandermonde potential Vydm, hence 
a = b = c = within numerical precision. This is shown 
in fig. 2. On the other hand, for (3 > (3c both coefficients 



B. Effective potential 

We are interested in the distribution of the eigenvalues 
Ai, A2 of the Wilson line L over the thermal ensemble of 
field configurations at each site. Here, we shall focus on 
the potential for the average of Ai and A2: 
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which for SU(2) is nothing but the absolute value of the 
Polyakov loop, p = yP". From the probability distribu- 
tion for p we define an effective potential via 



V eS (p) =- log P(p) 



(7) 



We normalize the probability density P{p) over the inter- 
val [0,1] which fixes the constant in the potential. This 
also factors out the volume dependence. The partition 
function can thus be written as 
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In our definition, therefore, V e s(p) is dimensionless as it 
absorbs the volume a 3 of a lattice cell and, implicitly, a 
factor of \/T (because we employ a 3d theory). 

Below, we shall show that in the weak-coupling regime 
(large (3, small m 2 ) a non-analytic contribution ~ Vi 2 to 
the effective potential arises dynamically. It is distinct 
from the Vandermonde potential VVdm = —\ log(l — £ 2 ) 
generated by the SU(2) group measure, and from the 



'bare" potential 



2 i 2 which is included in the ac- 



tion (5). We find that for a broad range of the couplings 
(3 and m 2 the potential has the form 



Vm(p) = -ilog(l-p 2 ) 



a — bp + cp 2 



(9) 



Note that the term proportional to p = V £ 2 of course 
does not break the Z(2) symmetry explicitly, and is not 
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FIG. 2: Effective potential for p = VP at m 2 = for various 
values of f3 fitted with the ansatz (9). For f3 < f3c (top) one 
sees only the group integration measure. Above (3c sponta- 
neous breaking of Z(2) is evident. 

a and b from eq. (9) turn non-zero. The quadratic coef- 
ficient c does not appear since m 2 — 0. The coefficient b 
appears to arise purely from the dynamics of fluctuations. 
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Its behavior is shown in fig. 3. 
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FIG. 3: Coefficients of the effective potential V(p) for m 2 — 
as a function of /3. Below /3c there is no deviation from the 
Vandermonde potential. Above /3c one sees a non-analytic 
contribution. (A discussion of the fit-curve is found in Ap- 
pendix A.) 



The ansatz (9) also works for m 2 ^ when j3 is not too 
large. We have confirmed this for a broad range of m 2 
for several fixed values of (3. Explicit results for j3 = 2.0 
are shown in fig. 4. The coefficient c of the quadratic 
term appears to coincide with the mass from the bare 
potential, c w — m 2 , within our present numerical pre- 
cision. However, we again observe that above the phase 
transition point the dynamics generates a non-analytic 
contribution to the potential for the Polyakov loop pro- 
portional to p = \J1? . For illustration, we show the m 2 
dependence of a, b, c for j3 = 1.0 and (3 — 2.0 in fig. 5. A 
detailed discussion of the dependence of a, b, c on (3 and 
to 2 is given in appendix A. The main point here is that 
the Vandermonde contribution to the effective potential 
does not depend on f3 and m 2 , and that the effective po- 
tential obtained from a histogram of p(x) at each lattice 
site is different from a Landau-Ginzburg type mean-field 
theory for the Polyakov loop. In the next section we shall 
see that when the field configurations are "cooled" to re- 
move short wavelength fluctuations, that in fact one does 
obtain a potential that resembles mean field theory, but 
with a coefficient multiplying the Vandermonde potential 
which depends on f3 and m 2 . 



C. Effective potential for block spins 

To obtain a potential for the long wavelength modes we 
average the Polyakov loop field over small cubes of side- 
length k before histogramming. This eliminates the short 
wavelength spatial field modes. We calculate blockspin 
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FIG. 4: Effective potential for p at (3 = 2.0 for various values 
of m 2 above and below the phase transition, fitted with the 
ansatz (9). 



averages defined as 
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(10) 



n= (0, 0, 0) . . . (jfe, k, k) 



We have investigated the cases k — 2,3,4. Blockspins 
were measured on 2500 independent configurations for 
each combination of (3, m 2 . As one expects from the cen- 
tral limit theorem, with increasing k the potential be- 
comes more symmetric and peaked around the actual 
expectation value (see fig. 7). In what follows, we take 
the configurations for k — 3 as a good approximation for 
the long distance sector. To see this consider fig. 6. At 
k = 3 the minimum of the fitted effective potential from 
eq. (11) differs at most by ~ 0.03 from the numerical 
result for the Polyakov loop expectation value (£}. This 
maximal deviation occurs exactly at the phase transition 
point. Here, within our numerical precision, k — 4 does 
not do significantly better. Away from the phase transi- 
tion, k = 3 differs less than ~ 0.01 from the numerical 
value of {£). k — 4 does slightly better but reduces our 
statistics significantly. All results presented below ap- 
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FIG. 5: Coefficients of V{p) for = 1.0 (top) and /3 = 2.0 
(bottom) as functions of m 2 . The non-analytic contribution 
sets in at the phase transition. 



pear to be stable when going from k — 3 to k = 4 (we 
will discuss an explicit example below). We refrain from 
discussing k — 2 in detail here as it appears that con- 
tributions from the short range fluctuation have not yet 
been completely eliminated. 

The ansatz (9) is no longer applicable for the long 
wavelength modes. The blockspin averaging appears to 
suppress the non-analytic term in the potential in most 
parts of the phase diagram (some possible exceptions are 
discussed below). It appears that over a broad range of 
f3 and m 2 the potential can be described fairly well by 
a form analogous to mean field models, with quadratic 
and quartic terms [18]. However, for a good fit to the 
extracted potential one needs to include an additional 
parameter do into the fit function which varies with /3 
and m 2 , and which multiplies the Vandermonde poten- 
tial. So, finally, the form which we use to model our data 
is 



V(p) = -d ±log(l 



di + d 2 p + d A p 



(11) 



FIG. 6: The difference of Polyakov loop expectation value 
and minimum of effective potential after "cooling" for k — 
2, 3, 4 at m 2 = 0.0 as a function of j3 (top) and at j3 = 2.0 
(bottom) as a function of m 2 . The peaks correspond to the 
phase transition point. 
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FIG. 7: Blockspin averages for p = 1.5, m 2 = 0.0. As k 
increases the potential peaks more sharply about the expec- 
tation value. 



Here, a linear term is not included. We also point out 
that the quartic term arises from the dynamics of fluc- 



tuations as such a contribution is not included in the 
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"bare" action (5). Explicit results for the numerical po- 
tential and for the fit via eq. (11), for different values 
of m 2 in the confined and deconfined phases and fixed 
/3 = 2.0 / 1.0, are shown in figs. 8,9. Below the phase 
transition one can set do — ^4 — and fit the poten- 
tial with a simple quadratic form. This is not surprising 
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FIG. 8: The potential for blocks of side length k = 3 fit 
with ansatz (11) for /3 = 1.0 in confined (top) and deconfined 
(bottom) phases. 

since for (£) = 0.0 the potential is essentially parabolic 
and there is little sensitivity to higher powers of I. 

We show how the fit parameters from eq. (11) depend 
on j3 and m? . Figs. 10 and 12 correspond to /3 = 1.0 with 
variable m 2 and m 2 = 0.0 with variable j3, respectively. 
One observes that right above the phase transition, the 
potential is a sum of quadratic and quartic terms while 
the Vandermonde contribution vanishes (in the figures, 
do appears to fluctuate somewhat around zero. We have 
checked, however, that the result is consistent with set- 
ting do = by hand). At higher m 2 or (3 respectively, the 
coefficient do increases gradually and saturates at about 
do rs 80 for large values of (3 or m 2 , while the quartic 
coefficient becomes negative. We have further checked 
that fixing do to its asymptotic value (do ~ 80) gives a 
less accurate fit and increases y 2 per degree of freedom in 
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FIG. 9: The potential for blocks of side length k = 3 fit 
with ansatz (11) for f3 = 2.0 in confined (top) and deconfined 
(bottom) phases. 
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FIG. 10: The coefficients in eq. (11) for /3 = 1.0. The Vander- 
monde term appears to gradually rise to its asymptotic value. 



the region closely above the phase transition roughly by 
a factor of two. This indicates that the gradual increase 
of do is a real dynamical effect and not just an artifact 
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p w 0.0 even slightly better than eq. (11), with a neg- 
ative coefficient d' Q in the region right above the phase 
transition. However, the improvement in \ 2 is below the 
percent level and the function (12) fails completely at 
large j3 or m 2 (by generating absurd global behavior). 

The function (12) may suggest that the suppression of 
the Vandermonde could be an artifact due to incomplete 
cooling of short-distance fluctuations. However, we have 
investigated the cases [/3 = 1.0/variable m 2 ] and [m 2 = 
0.0/variable f3] also for k = 4 and obtained similar results, 
up to an overall scaling factor for all coefficients in the 
potential 3 . We show the result for j3 = 1.0 in fig. 13. 
Compared to fig. 10 we only observe a slight suppression 
of the quartic coefficient right above the phase transition. 



FIG. 11: The coefficients of eq. (11) (3 = 2.0. The region with 
vanishing Vandermonde term appears to shrink when going 
from = 1.0 to = 2.0. Compare to fig. 10 
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FIG. 13: The coefficients of eq. (11) for /3 = 1.0 at k = 4. 
Compare to fig. 10. 



FIG. 12: A region with vanishing Vandermonde term is visible 
for m 2 = 0.0 and ,3 > Re- 
generated by a lack of sensitivity to the Vandermonde 
potential when (p) is much smaller than 1. Fig. 11 shows 
that the region right above the phase transition, where 
the Vandermonde vanishes, appears to shrink when going 
deeper into the weak coupling limit (larger (3). 

The question remains, how this suppression of do just 
above the phase transition, at moderately weak coupling, 
comes about. We have therefore attempted to model 
the potential right above the transition with a differ- 
ent function, assuming a fixed Vandermonde potential 
term, equal to the asymptotic value, but also allowing 
additional terms. We find that it is possible, within our 
numerical accuracy, to trade the suppression of the Van- 
dermonde for another term linear in p. The function 

V(p) = -d \ log(l - p 2 ) + di + d' p + d 2 p 2 + d 4 p 4 , 

with d = 80 . (12) 

reproduces the behavior of the effective potential around 



IV. SUMMARY AND DISCUSSION 

We have performed simulations of an effective theory 
of Wilson lines coupled to gauge fields in three dimen- 
sions which respects the center symmetry of the four- 
dimensional SU(2) Yang-Mills theory. After mapping the 
phase diagram, we have investigated the effective poten- 
tial for the average of the eigenvalues of the SU(2) Wilson 
line, which is equal to the absolute value of the Polyakov 
loop. We found that a form containing non-analytic con- 
tributions can describe the extracted potential. This non- 
analytic term was not present in the action, and therefore 
must arise from the dynamics. 



3 A rescaling of the coefficients in the potential is expected be- 
cause coarse-graining over k lattice sites effectively corresponds 
to a simulation with different lattice spacing. A quantitative 
comparison of the coefficients of the k = 3 potential to those of 
the k = 4 potential would also require a rescaling of (3. 
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We extracted a similar effective potential also for the 
long wavelength modes of the Polyakov loop and found 
that this can be described by a mean-field type poten- 
tial with quadratic and quartic terms plus an effective 
Vandermonde potential which depends on the couplings. 
Just above the phase boundary, in the deconfined phase, 
the effective Vandermonde potential contributes little. 
Deeper into the deconfined phase its coefficient increases 
and eventually appears to approach a constant. 

Our simulations may provide useful insight into the 
structure of mean-field type models for the deconfining 
phase transition. For example, so-called "Polyakov-NJL" 
models have recently been studied extensively. Such 
models attempt to describe QCD thermodynamics over a 
range of quark masses, from the pure-gauge limit to phys- 
ical QCD; they require an ansatz for the effective poten- 
tial for the Polyakov loop. For example, in early works on 
this subject [19] a quadratic potential for I has been used, 
plus a Vandermonde contribution (per lattice site) which 
is constant and temperature independent. Our results 
appear to indicate, however, that if a standard potential 
with terms ~ C 1 and ~ £ 4 (plus cubic Z(3) invariants 
for the case of three colors) is used, that a temperature 
dependent Vandermonde contribution should also be al- 
lowed for. 
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with (3c = log(l + V2) (see fig. 5). This is similar to the 
magnetization in the 2D Ising model [20] , which is given 
by 



M(/3) = 6(0 - (3 C ) (1 - (sinh(2/3J)- 4 )s; 



(A3) 



We find that (A2) works reasonably well for j3 < 3.0 (see 
the fit-curves in fig. 5 for specific examples). We have in- 
cluded the constant (3c in eq. (A2) because it corresponds 
to the critical point of the Ising model. Our model of 
course deconfines at a different value of (3. Nevertheless, 
an ansatz such as eq. (A3) implicitly assigns the number 
(3c a special meaning and we therefore include it into our 
ansatz also in order to "filter out" its effect. Isolating ftc 
in such a way simplifies the resulting dependence of the 
fit parameters on (3 greatly. 

The coefficients introduced in eq. (A2) act as follows: 
m 2 corresponds to a shift along the horizontal axis, b is 
a scaling factor and g is a modification of the coupling 
strength. The (3 dependence of these coefficients can be 
described by power laws 



m 2 {(3) 
b((3) 

m 



ml + m' 2 (3 v 
b' + b' '(3 v , 
.90 + g'o0 u , 



(A4) 
(A5) 
(A6) 



with 



= 2.2(1) , m 2 = -2.1(1 



b'o = 
9o = 



-2.5(4) , b' ' = 6.0(4), v 
0.038(1) , 5o = 0.017(l) , 



w = -1.2(1) ,(A7) 
= 0.90(4) , (A8) 
u = 2.8(1) . (A9) 



The (3 dependence of m 2 , b and g is shown in fig. 14. 



APPENDIX A: COEFFICIENTS IN THE 
EFFECTIVE POTENTIAL FOR ALL MODES 



In this appendix we discuss the behavior of the non- 
analytic term in the potential defined in eq. (9) as a func- 
tion of f3 and m 2 . 

The case m 2 = 0.0 is rather simple. Fig. 3 suggests an 
ansatz of the form 

609) 



b (f3 - pcY 0(P - Pc) ■ 



(Al) 



Indeed, we find that with (3c = 0.9, a good fit of b((3) is 
possible, resulting in b = 7.1(1) and r = 0.82(4). This 
fit corresponds to the solid line in fig. 3. 

For m 2 ^ 0.0 the situation is more involved. Motivated 
by the 2D Ising model we try the ansatz 

b((3,m 2 ) = b((3)8{m 2 +rh 2 {(3)) 



sinh g(j3) (m 1 + to (/?)) + j3 c 



(A2) 
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FIG. 14: f3 dependence of various coefficients parameterizing 
b((3, m 2 ) with their corresponding fit curves. g(/3) is scaled up 
by a factor of 30. 
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APPENDIX B: FIRST-ORDER TRANSITION IN 
EXTREME WEAK COUPLING LIMIT 

At very large (3 the transition becomes first order. To 
see this, consider figs. 15 and 16. For f3 = 5.0 the effec- 
tive potential develops two distinct minima in the vicin- 
ity of the phase transition point. Correspondingly, the 
Polyakov loop expectation value appears to be discon- 
tinuous. This behavior is in sharp contrast to the case 



(3 = 2.0 (shown in fig. 4), where one can see a single 
minimum moving continuously with m 2 (within the res- 
olution) from p w to p w 1 . 

The first order transition indicates that our effective 
theory cannot describe 4D Yang Mills when (3 is too large. 
It is possible that this behavior is cured by adding higher 
powers of the Polyakov loop to the action (5). 
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FIG. 15: Effective potential for f3 = 5.0 at values of m 2 
slightly above and slightly below the phase transition. A first 
order phase transition is apparent since one observes two dis- 
tinct minima. 



FIG. 16: Polyakov Loop expectation value for = 2.0 and 
f3 — 5.0 measured on N s = 24. The transition becomes very 
sharp for large j3. 
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